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ABSTRACT 

We present stellar velocity dispersion profiles for seven Milky Way dwarf spheroidal (dSph) satellite 
galaxies. We have measured 8394 line-of-sight velocities (±2.5 km s _1 ) for 6804 stars from high- 
resolution spectra obtained at the Magellan and MMT telescopes. We combine these new data with 
previously published velocities to obtain the largest available kinematic samples, which include more 
than 5500 dSph members. All the measured dSphs have stellar velocity dispersion of order 10 km 
s _1 that remains approximately constant with distance from the dSph center, out to and in some 
cases beyond the radius at which the mean surface brightness falls to the background level. Assuming 
dSphs reside within dark matter halos characterized by the NFW density profile, we obtain reasonable 
fits to the empirical velocity dispersion profiles. These fits imply that, among the seven dSphs, 
M V i r ~ 1O 8 ~ 9 M . The mass enclosed at a radius of 600 pc, the region common to all data sets, 
ranges from (2 - 7) x 1O 7 M . 

Subject headings: galaxies: dwarf — galaxies: kinematics and dynamics — (galaxies:) Local Group 
— (cosmology:) dark matter — techniques: radial velocities 



1. INTRODUCTION 

Spanning the absolute magnitude range —13 < 
My < -4 (|Mateol [1991 iBelokurov et al.ll2007t ). dwarf 
spheroidal (dSph) galaxies represent the extreme end 
of the galaxy luminosity function. The Milky Way's 
(MW's) ~ 15 known dSph satellites have managed thus 
far to avoid and/or to survive interactions that consumed 
many of their siblings during hierarchical formation of 
the MW halo. Recent veloci ty measurements for sev- 
eral hundred stars per dSph dKlevna et al.1 l2002t I2003J . 
2001 iMunoz et al.l I2005L 120061 : iWalker et all 12006^ 



Simon &: Gehal l2007f ) demonstrate that dSph velocity 



dispersions remain approximately flat with radius. Anal- 
yses that recover M(r) from the Jeans equation conclude 
that dSphs poss ess extended and dominant DM halos 
(e.g.. lLokasll2002D . unless stellar velocity distributions are 
highly anisotropic or ongoing tidal d isruption inva lidates 
the assumption of equilibrium (e.g.. lKroupalll997l ). 

In this letter we present velocity dispersion profiles for 
seven dwarf satellites of the Milky Way — Carina, Draco, 
Fornax, Leo I, Leo II, Sculptor, and Sextans. Profiles 
are calculated from new velocity data we have obtained 
using the Magellan and MMT telescopes (| Walker et al.1 
I2007t iMateo et all [20071 ). After combining with previ- 
ously published data, samples range in size from ~ 200 
(Leo II) to ~ 2000 (Fornax) member stars per dSph, for 
a total of 5544 probable members in the seven dSphs. 
The new samples more than double the amount of ex- 
isting dSph kinematic data and allow us to measure ve- 
locity dispersion profiles with unprecedented precision. 

"THIS PAPER INCLUDES DATA GATHERED WITH THE 6.5- 
M MAGELLAN TELESCOPES AT LAS CAMPANAS OBSER- 
VATORY, CHILE, AND WITH THE 6.5-M MMT TELESCOPE. 
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All the measured dSphs exhibit approximately flat ve- 
locity dispersion profiles that can be fit reasonably well 
under the assumption that dSphs are equilibrium sys- 
tems embedded within dark matter halos. We show that 
if the dark mat ter halos are charact erized by the NFW 
density profile ([Navarro et al.l [1996) that results from 
N-body simulations of cold dark matter, the kinematic 
data suggest dSphs have M vir ~ 1O 8_9 M and masses 
of (2 - 7) x 10 7 Af Q inside a radius of 600 pc. 

2. OBSERVATIONS & DATA 

Over seven observing runs between 2004 March and 
2007 January, we used the Michigan/MIKE Fiber Sys- 
tem (MMFS) at the Magellan/Clay 6.5m telescope to 
obtain 7383 high-resolution echelle spectra of 5793 in- 
dividual red giant candidates in th e dSph galaxies Ca - 
rina, Fornax, Sculptor and Sextans. I Walker et al.l (2007) 
describe the uniform data reduction procedure in de- 
tail. MMFS spectra sample the Mg-triplet (MgT) re- 
gion, spanning 5140 — 5180 A with effective resolution 
- 0.1 A/pix (R ~ 20000). We measure both stellar veloc- 
ity and the pseudo-equivalent width of the magnesium- 
triplet absorption feature. The latter quantity correlates 
with surface gravity and helps determine dSph member- 
ship. Comparisons of repeat measurements for more than 
1000 stars indicate median measurement errors of ~ 2.0 
km s _1 and ~ 0.06 A, respectively. 

Additionally we obtained 1183 spectra from 1011 red 
giant candidates in the dSphs Leo I, Leo II, and Draco 
using the multi-fiber Hectochcllc spectrograph at the 
MMT 6.5m telescope during three obser ying runs in 
Marc h/April of 2005, 2006 and 2007 (see IMateo et all 
|2007| ). Hectochelle spectra sample the MgT region over 
5150-5300 A with effective resolution 0.01 A/pix (R ~ 
25000). Repeat measurements for more than 100 stars 
with Hectochelle indicate velocity errors of ~ 2.6 km s~ x . 

To these new, homogeneous data sets we add 2239 
velocities previously published for red giant candi- 
dates in the observed dSphs (Carina: IMateo et al.l 
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1991 IMufioz et all [20061: D raco: lArmandroff et al" 



199a IKlevna et all l2002t Forn ax: IMateo et al 



199li IWalker et alj l2006al: Leo I: IMateo et all 1199 



TABLE 1 

Summary of Velocity Samples and NFW Parameters 



Koch et al.l l2007U Le o ILlVogt et a ij |1995l : IKoch et all 
2007a|; Sculptor: EVestfall et alTl2006l : Sextans: 
Hargreaves et al.lfl99l: IKlevna et al.ll2004l ). We correct 



for zero-point offsets and combine measurements of com- 
mon stars using the weighted (by measurement errors) 
mean velocity . We d o not correct for the slopes noted by 
IWalker et al.l (|2007f) when comparing MgT velocities to 
those measured from the calcium triplet; because the rel- 
evant combined samples are dominated by MMFS data, 
doing so would have no significant impact on the velocity 
dispersion profiles. 

We evaluate the probability, Pdsphi , that the i th star is 
a dSph member using up to three pieces of information — 
velocity VI, magnesium strength W, (available only for 
MMFS samples), and position — a priori, stars projected 
farther from the dSph center are less likely to be mem- 
bers. We assume 1) the joint distribution of V and W 
for the members of a given dSph is a bivariate Gaus- 
sian; 2) interlopers have magnesium strengths following 
a univariate Gaussian distribution; and 3) interlopers 
have a non-Gaussian velocity distribution that we esti- 
mate numerically f rom the Besangon Milky Way model 
(jRobin et a l. 2003) . We use an iterative, expectation- 
maximization (EM) algorithm ( (Sen et aTl |2007| ) to re- 
cover maximum-likelihood estimates of the Gaussian 
means and variances as well as the individual member- 
ship probabilities (see Walker et al. in prep, for full de- 
tails of the algorithm). Black markers in FigureQ]idcntify 

the most likely interlopers (Pdsphi < 0.01), while colored 

markers signify larger Pdsphi ( see caption to Figure Q]) . 
In subsequent calculations we weigh each data point by 

Pdsphi ' 

Middle sub-panels in Figure Q] plot velocity against an- 
gular distance from the dSph center. Bottom sub-panels 
indicate the overall velocity distribution; red lines give 
the expected distributions of interloper velocities, per 
the Besangon model. For dSphs observed with a high 
degree of contamination (Carina, Sextans), we find ex- 
cellent agreement between the predicted and observed 
velocity distributions of foreground stars. Table [T] lists 
for each dSph the numbers of newly observed stars, total 
stars after combining with published data, and member 

stars Ndsph = Yn=i Pdsph, ■ 

3. VELOCITY DISPERSION AND NFW PROFILES 

We estimate line-of-sight velocity dispersion profiles af- 
ter dividing each sample into ~ */ Nd sp h bins accord- 
ing to projected distance from the dSph center. For a 
given sample we define bins (circular annuli) such that all 
contain approximately equal numbers of dSph members. 
Thus the number of stars, including interlopers, in each 
bin may vary but for all bins, Pdsphi ~ \/N dsph . 

We use a Gaussian maximum-likelihood method (see 
IWalker et al.ll20"06af ) to estimate the velocity dispersion 
within each bin. 

Left-hand panels Figure [5] display the resulting veloc- 
ity dispersion profiles, which generally are flat. The 
outer profile of Draco shows no evidence for a rapidly 
falling dispersio n , con trary to evidence presented by 
Wilkinson et alj ((2004), but consistent with the result 
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of IMufioz et all (|2005| ) 5 . In fact the outer profiles of 
Draco, Carina and perhaps Sculptor show gently ris- 
ing dispersions. While it is likely that at least in Ca- 
rina thi s behavior is assoc ia ted with the onset o f tida l 
effects ()Munoz et al.ll2006D . iMcConnachie et al.1 (|2006h 
point out that the tendency of some dSphs to have sys- 
tematically smaller velocity dispersions near their cen- 
ters is perhaps the result of distinct and poorly mixed 
stellar populations dTolstov et all 120041 : iBattaglia et al.l 
2006; Ibata et al. 2006). Either explanation complicates 
a thorough kinematic analysis; in the present, simplified 
analysis we assume all stars belong to a single population 
in virial equilibrium. 

Dashed lines in Figure[2]are velocity disper sion profiles 
calculated for single-component King models (King 1962) 
conventionally used to characterize dSph surface bright- 
ness profiles. The adop t ed Ki ng models are those fit by 
llrwin fe Hatzidimitrioul (|1995l hereafter IH95) and nor- 
malized to match the central velocity dispersion. Their 
systematic failure to predict flat profiles removes any 
doubt that (Newtonian) mass- follows- light, equilibrium 
models provide a poor description of dSph kinematics. 

In forthcoming papers we examine evidence for tidal 
effects and explore the range of dark matter density pro- 
files that are consistent with the dSph velocity data. 
Here we simply show that dark matter halos having 
the NFW density profile can fit the data reasonably 
well. Under the assumptions of 1) spherical symme- 
try, 2) dynamic equilibrium, 3) radially constant veloc- 
ity anisotropy (3 = 1 — cr| /of , and that 4) dSph surface 
brightness declines exponentially with radius, with cen- 
tral value and scale length adopted fr om IH95, we use 
the J eans equation (e.g., Eq. A15 of iMamon fc Lokasl 
2005) to calculate projected velocity dispersion profiles 
expected for NFW halos over a range of plausible virial 
masses. Assuming the rela tionship between M V j r and 
halo concentration found bvlJing fc Sutol ([2000L see also 
IMamon fc Lokasl [20051: IKoch et al.ll2d07bh . we fit to the 
empirical velocity dispersion profiles using just two free 
parameters, M V i r and (3. Solid lines in the left-hand pan- 
els of Figurc[2]indicate the best-fitting NFW models, and 
Table[T]lists the associated parameters, including M rmaic , 
the enclosed (dark plus luminous) mass at the outermost 
profile point. The right-hand panels of Figure [2] display 
mass-density, mass and M/L profiles corresponding to 
t he best- fitting NFW m odels. 

IStrigarT et all (|2007f l show that dSph masses within 
r < 600 pc (Mgoo) are constrained robustly by kinematic 



5 W e h ave not included the unpublished data of Wilkinso n et al.l 
pOOl ) or [Munoz et al. ( 2005) in our calculations of the velocity 
dispersion profiles of Draco. 
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Fig. 1. — Mg pseudo-equivalent width, projected distance, and membership fraction vs. velocity for new dSph samples obtained with 
Magellan/MMFS and MMT/Hectochelle (Mg strength is unavailable for Hcctochelle observations). Marker color indicates the probability 
that the star is a dSph member. Black (red; magenta; green; cyan; blue) markers signify Pdsph < 0-01 (Pdsph > 0.01; > 0.50; > 0.68; > 
0.95; > 0.99). Pairs of dotted lines enclose stars that would pass conventional membership tests based on a 3<r velocity threshold. Red lines 
in the bottom sub-panels indicate the interloper distribution predicted by the Besangon Milky Way model. 



data, irrespective of the assumed form of the density 
profile. This radius is convenient for practical reasons 
in our analysis as it covers the largest region common to 
all the data sets and is sufficiently small that one expects 
tidal effects inside this radius to be negligible (Re ad et al.l 
l2006f ). Integrating the density profiles for the stellar (for 
which we assume M/L = 1) and dark components, we 
find that Mqqq ranges from (2 — 7) x 1O 7 M among the 



seven dSphs (Table [T]). This narrow range is consistent 
with the notion that these dSphs, while exhibiting order- 
of-magnitude va riability in luminos ity, occupy dark halos 
of similar mass (jMateo et al.lll993l ). 

We thank the staff at the Las Campanas and MMT 
Observatories. This work is supported by NSF Grants 
AST 05-07453, AST 02-06081, and AST 94-13847. 
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